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Abstract 

Using data from Fermilab fixed-target experiment E791, we have measured particle- 
antiparticle production asymmetries for A , H - , and Q~ hyperons in 7r~ - nucleon 
interactions at 500 GeV/c. The asymmetries are measured as functions of Feynman- 
x (x F ) and p% over the ranges —0.12 < x F < 0.12 and < p% < 4(GeV/c) 2 . We 
find substantial asymmetries, even at x F = 0. We also observe leading-particle-type 
asymmetries which qualitatively agree with theoretical predictions. 



Strange particle production is an important tool for studying how non-perturbative 
QCD affects light quark production and hadronization. One process that in- 
volves both production and hadronization is the leading particle effect. This 
effect is manifest as an enhancement in the production rate of particles which 
have one or more valence quarks in common with an initial state hadron 
compared to that of their antiparticles which have fewer valence quarks in 
common. This enhancement increases as the momentum of the produced par- 
ticle increases, in the direction of the initial hadron with which the produced 
particle shares valence quarks. This process has been extensively studied in 
charm production in recent years from both experimental [1,2] and theoret- 
ical [3] points of view. The same type of leading particle effects have been 
seen in light hadron production [4] and are expected [5]. Other effects, like 
the associated production of a kaon and a hyperon, can also contribute to an 
asymmetry in hyperon- antihyperon production [6]. 

As a byproduct of our charm program in Fermilab experiment E791, we col- 
lected a large sample of A — A , 2~ — S + , and Q~ — Q + hyperons which 
we have used to measure the particle-antiparticle production asymmetries re- 
ported here. For hyperons, the range of scaled longitudinal momentum x F 
(2p L /y/s) covered by our experiment was —0.12 < x F < 0.12. Given E791's 



negative pion beam and nucleon target, we expect recombination effects [5] 
to produce different asymmetries in the beam and target fragmentation re- 
gions as the content of valence quarks in the beam and target hadrons differ. 
A growing asymmetry is expected in A — A production as \x F \ increases 
in the negative direction towards the target fragmentation region. A smaller 
asymmetry (or none at all) is expected for x F > 0, towards the beam frag- 
mentation region, because A and A each share one valence quark with the 
incident pion. For H~ — S + production, a growing asymmetry with \x F \ is 
expected in both regions since the S~ shares one valence quark with the 7r~ 
as well as with the target particles (p and n) whereas the S + shares none. 
No leading particle effects are expected for Q~ or Q + as they have no valence 
quarks in common with either the target or the beam. 

Measurements of A — A [7], 2~ — S + , and Q,~ — f2 + [8] production asymme- 
tries in 7r~Cu interactions at 230 GeV/c have been reported in the literature. 
Additional measurements of the A ° — A asymmetry can be found in the ex- 
periments listed in [9], but in general light hadron production asymmetries in 
n~N interactions have not been studied extensively. 

Experiment E791 recorded data from 500 GeV/c 7r~ interactions in five thin 
foils (one platinum and four diamond) separated by gaps of 1.38 to 1.39 cm. 
Each foil had a thickness of approximately 0.4% of a pion interaction length 
(0.5 mm for the upstream platinum target, and 1.6 mm for each of the carbon 
targets). The E791 spectrometer [10] in the Fermilab Tagged Photon Labo- 
ratory was a large-acceptance, two-magnet spectrometer augmented by eight 
planes of multiwire proportional chambers (MWPC) and six planes of silicon 
microstrip detectors (SMD) for beam tracking. The magnets provided a total 
transverse momentum kick of 512 MeV/c. Downstream of the target there 
were 17 planes of SMD's for track and vertex reconstruction, 35 drift chamber 
planes, two MWPC's, two multicell threshold Cerenkov counters, electromag- 
netic and hadronic calorimeters (with apertures about 70 by 140 mr), and a 
muon detector. An important element of the experiment was its extremely fast 
data acquisition system [11] which was combined with a very open transverse- 
energy trigger to record a data sample of 2 x 10 10 events. The trigger required 
that the total "transverse energy" (i.e., sum of the products of energy observed 
times the tangent of the angle from the target to each calorimeter segment) 
be at least 3 GeV. 

For this analysis we use only interactions in the isoscalar carbon targets so 
that our results truly represent a "nucleon" , that is, the average of neutrons 
and protons. Most A°'s decay before entering the drift chamber region (150 
cm downstream of the targets) but downstream of the end of the silicon vertex 
detectors (50 cm from the targets), while some S^'s and Q _ 's decay in the 
silicon region. 



Throughout this paper, references to a particle should be taken to include its 
antiparticle unless explicitly stated otherwise. A°'s were reconstructed using 
the p7r~ decay mode. Proton and n~ tracks were required to have a distance 
of closest approach less than 0.7 cm at the decay vertex and to have an in- 
variant mass between 1.101 and 1.127 GeV/c 2 . In addition, the ratio of the 
momentum of the proton to that of the pion was required to be larger than 
2.5. The reconstructed A decay vertex formed by the two tracks was required 
to be downstream of the last target but upstream of the first magnet. For the 
A production study, we removed A°'s coming from S~ decay by requiring 
that the A ° candidates have an impact parameter with respect to the primary 
vertex of less than 0.3 cm if decaying within the first 20 cm downstream of 
the target, and less than 0.4 cm otherwise. After these cuts, the remaining 
2~ contamination was ~ 1.5%, having negligible effect on the A — A asym- 
metry. The reconstructed mass distributions for this sample, for each interval 
of x F and p 2 , were fit using a binned maximum likelihood method with a 
Gaussian signal plus a linear background. All bin widths were much larger 
than the experimental resolution of the variable binned. In the fit, the central 
reconstructed mass values and mass resolutions were fixed to values obtained 
from Monte Carlo simulation. Contamination from misidentified K® — > 7t + tt~ 
decays provided a flat background in these distributions in the range shown, 
and produced a negligible effect on the fit numbers of A°'s. The reconstructed 
mass distributions for the entire sample are shown in Figs. 1(a) and (b). For 
the A —A analysis we use data from approximately 7% of the overall sample 
recorded for the experiment. The total signal, taken as the sum of background 
subtracted signal in each bin, was 2 587870 ± 1780 A°'s and 1690 030 ± 
1500 A°'s. 



2~ were selected via the decay mode H~ — » A°7r~ and Q~ via the decay mode 
Q~ — ► A°K~. Beginning with a A candidate, we added a third, distinct 
track as a possible pion or kaon daughter. All three tracks were required to 
have hits in the drift chamber region only. For these samples, we removed the 
requirement on the A ° impact parameter. The invariant mass of the candidate 
hyperon (calculated from the known mass and measured momentum of the 
A , as determined by its two decay tracks, together with the third track) 
was required to be between 1.290 and 1.350 GeV/c 2 for the H - and between 
1.642 and 1.702 GeV/c 2 for the fi~ . The charged hyperon track, entirely in 
the SMD region, was required to be reconstructed. Its direction had to match 
that determined by its daughter tracks in the drift chamber region to within 1 
mrad and it was required to have a distance of closest approach to the primary 
vertex of less than 100 /xm. The S~ or Q~ vertex (reconstructed from the A 
and daughter tt~ or K~) was required to be downstream of the last SMD plane 
and upstream of the A vertex. The former requirement allows us to track the 
hyperon in the silicon. For the Q~ sample, the third track was required to have 
a kaon signature in the Cerenkov counters and momentum in the range 6-40 



GeV/c. For this momentum range, the Cerenkov kaon identification efficiency 
was about 85% and the pion misidentification rate was about 5%. 

As with the A sample, the S - — S + and Q~ — Q + invariant mass plots were 
fit to a Gaussian signal plus linear background for each interval of x F and p\. 
The total reconstructed mass distributions are shown in Figs. 1(c) through (f). 
With the final selection criteria and the full E791 data set, we found 996 180 
± 1200 E~, 706 620 ± 1020 5+ 8 750 ± 110 Q", and 7460 ± 100 fi+ after 
background subtraction. Again, these numbers and their errors come from the 
sum of signals in all bins. We checked that the S~ contamination for the fi~'s, 
after all cuts, was negligible. 



For each x F and p\ bin and for each hyperon Y, we defined an asymmetry 
parameter A as 

A _ N v - N- rv e v 

A v _ v = 



N Y + N-r Y ' e-' 

where N Y (N-) is the number of hyperons (antihyperons) produced in the 
bin, and e Y (e— ), the product of the geometrical acceptance and reconstruction 
efficiency for each hyperon (antihyperon). Values for the iV's were obtained 
from the individual fits to the mass plots for events selected to lie within each 
x F and p\ range. 

Selection criteria for the particle and antiparticle samples were identical. How- 
ever, geometrical acceptances and reconstruction efficiencies were not neces- 
sarily the same, mostly due to an inefficient region in the drift chambers 
produced by the negative pion beam. To evaluate this effect, a large sample of 
Monte Carlo (MC) events was created using the Pythia/Jetset event gen- 
erators [12]. These were projected through a detailed simulation of the E791 
spectrometer to simulate "data" in digitized format which was then processed 
through the same computer reconstruction code as that used for data from the 
detector. Candidate events were then subjected to the same selection criteria 
as that used for data. To account for correlations between x F and p T , efficien- 
cies were determined in bins of the two parameters. The ratios of efficiencies 
e Y /e— (x F ,p F ) evaluated from the MC samples are shown for each hyperon in 
Figs. 2(a)-(c). For the x F range of our data, there are only minor differences 
between particle and antiparticle efficiencies. This result was expected, as hy- 
peron decays with tracks aimed at the inefficient part of the drift chambers 
were rare. As a further check, we verified that our results were not significantly 
affected by symmetrically eliminating these events. 

The asymmetries obtained using Eq. 1 are shown in Figs. 3 and 4 for the 
hyperons A , E~, and H", as a function of x F and p\, respectively. For each 



hyperon the asymmetry is presented as functions of x F for different intervals 
of p%, and also as functions of p\ for different intervals of x F . The asymmetries 
show a substantial dependence on x F and some dependence on p\. There is 
evidence for a correlation between x F and p\ 



'T- 



The asymmetry A(x F ) integrated over p^,, and the asymmetry A(jj%) integrated 
over x F , are presented in Fig. 5. This figure shows asymmetries for all three 
types of hyperons. The results are also listed in Table 1 along with statistical 
and systematic errors. 

We looked for systematic effects from the following sources: 

• Event selection criteria; 

• The minimum transverse energy in the calorimeters required in the event 
trigger; 

• Uncertainties in calculating relative efficiencies for particle and antiparticle; 

• Effect of the 2.5% K~ contamination in the beam. 

The first two effects were found to be negligible. A significant error came 
from the spectrometer efficiencies, and these uncertainties are included as 
systematic errors in Table 1. 

The effect of the K~ contamination in the beam was difficult to estimate as no 
data on hyperon asymmetries in K ~ production in this x F range existed. How- 
ever, even for equal n~ and K~ production of A°'s (as opposed to A 's) and 
100% asymmetry for kaons (i.e., no A production), only a 1.5% change would 
occur in the final result. As for the 7i~ (ud) beam, no rise in the asymmetry 
in the positive x F direction was expected since both A°(uds) and A (uds) 
shared one valence quark with the K~(us) beam. Similarly, any rise in the 
^(dss) — S + (dss) asymmetry should have been approximately the same as 
for a 7r~ beam since, in both cases, the S~ hyperon shared one valence quark 
with the beam particle. Since the f2~(sss) was a leading particle for a K~ 
beam, in this case some asymmetry was expected in the forward x F region. 
We noted the low K~ content of the beam and the low x F of our data. In any 
event, no effect was evident. Thus we did not correct our measured values of 
asymmetry for possible effects of kaon contamination in the beam. 

The asymmetry curves for A and for H~ cross over at x F ~ 0.02. At larger 
positive x F the H~ asymmetry becomes considerably larger than that for A . 
The fact that both beam and target particles share one valence quark with 
the S~ and none with the S + , would predict that the S~ asymmetries rise 
with increasing \x F \ in both positive and negative x F regions, as observed 
in our data. In contrast, the A shares two valence quarks with the target 
and only one with the beam particles, while the A only shares one with the 
beam. Consequently the A — A asymmetry rises sharply for x F < and 



not for x F > 0. The Q~ — Q + asymmetry appears to be constant. Notice that 
neither the Q~ nor the Q + shares valence quarks with either the beam or 
target particles. 

The observed positive asymmetry for Q could arise, in part, from the associ- 
ated production of kaons. This could also be the origin of the observed positive 
asymmetry near x F = 0.0 for S's and A's. Part of the associated production 
enhancement of baryons, as opposed to antibaryons, may come from the higher 
energy thresholds for the production of antibaryons. For production of hyper- 
ons, the conservation of strangeness requires only the associated production of 
one or more kaons. For production of an antihyperon, baryon number as well 
as strangeness must be conserved, requiring the associated production of at 
least two additional baryons, thus raising the energy thresholds and favoring 
particle over antiparticle production [6]. 

The behavior of the three asymmetries shown in Fig. 5 gives evidence for the 
leading particle effect. In the backward region (x F < 0), a larger asymmetry is 
observed when there is a larger difference in the number of valence quarks in 
common with the target, for hyperon and antihyperon. Evidence for a similar 
effect in the production of D ± mesons in the forward region (x F > 0) was 
presented by the E791 collaboration [1] and others [13]. The E791 collabora- 
tion has also studied asymmetries in the production of Df [2] and A^ [14]. 
However, leading particle effects for charmed hadrons typically occur at larger 
values of x F than they do for strange hadrons. 

The Pythia/Jetset [12] model describes only some features of our results, 
and those only qualitatively, as can be seen in Fig. 5. This model predicts 
small values of asymmetry for x F ~ 0; this is in contrast with our results 
which range from 0.08 to 0.18 in this region. Our data show that, even in the 
central region, the asymmetries are not zero, and suggest that leading particle 
effects play an even larger role than expected as \x F \ increases. 

In summary we report the most precise, systematic study to date of the pro- 
duction asymmetry for A, S, and Q hyperons in a single experiment. The range 
of x F covered, —0.12 < x F < 0.12, allows the study of asymmetries in regions 
close to x F = for the first time in a fixed target experiment. Some evidence 
for possible correlations between x F and p% is observed (see Figs. 3 and 4). 
Our results for particle-antiparticle asymmetries are consistent with the ex- 
perimental results obtained by previous experiments [7,8] (see Table 2) but 
with smaller uncertainties. Our results can be described qualitatively in terms 
of the energy thresholds for the production of hyperons and antihyperons to- 
gether with their associated particles and a model in which the recombination 
of valence and sea quarks in the beam and target particles contributes to the 
hyperon and antihyperon production in an asymmetrical manner [5]. 
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Fig. 1. Effective mass distributions for decay products of hyperons in events selected 
for this study and the corresponding signals with background subtracted. Plots 
correspond to (a) A — ► p7r~; (b) A — ► p7r + ; (c) H~ — > A°7r~; (d) H + — ► A 7r + ; 
(e) Q~ — ► A i^ _ (f) f2 + — ► A X + . The numbers and errors in the top corner of 
each plot come from the sum of the number of background subtracted events in the 
individual {x F ,p\) bin fits. 
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asymmetries as a function of pj, in different x F regions. 
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Fig. 5. Comparison among the asymmetries in the production of the three hyperons 
for the E791 data and Pythia/ Jetset predictions. The error bars of the E791 data 
include the statistical and systematic uncertainties 
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Table 1 

Final asymmetries of A, E, and 0, showing the statistical plus systematical errors. 
The statistical errors given include those due to the number of observed events and 
the number of MC events. Typically, the error is dominated by the MC uncertainty. 

x F Region A AO _^o = ^^ A s -_ s+ = §-^ A n -_ n+ = g-^ 

-0.12 - -0.09 0.392 ± 0.007 ± 0.001 0.286 ± 0.014 ± 0.001 0.047 ± 0.096 ± 0.028 

-0.09 - -0.06 0.315 ± 0.005 ± 0.001 0.220 ± 0.010 ± 0.001 0.128 ± 0.046 ± 0.005 

-0.06 - -0.03 0.240 ± 0.004 ± 0.001 0.182 ± 0.007 ± 0.001 0.100 ± 0.030 ± 0.001 

-0.03 - 0.00 0.179 ± 0.002 ± 0.001 0.154 ± 0.007 ± 0.001 0.099 ± 0.027 ± 0.002 

0.00 - 0.03 0.136 ± 0.002 ± 0.001 0.139 ± 0.009 ± 0.001 0.106 ± 0.028 ± 0.001 

0.03 - 0.06 0.111 ± 0.005 ± 0.001 0.158 ± 0.012 ± 0.001 0.085 ± 0.038 ± 0.005 

0.06 - 0.09 0.094 ± 0.008 ± 0.001 0.196 ± 0.017 ± 0.001 0.075 ± 0.058 ± 0.011 



0.09 - 0.12 


0.086 ± 0.012 ± 0.003 


0.250 ± 0.025 ± 0.001 


0.129 ± 0.106 ± 0.018 


p% Region 


a A°-A° 
A A°-A° A o+A u 


a - E--S+ 


A Q--C1+ 

An--n+ - n - +n + 


■ - -- H-+S+ 



0.0 - 0.5 0.206 ± 0.002 ± 0.001 0.162 ± 0.007 ± 0.001 0.094 ± 0.034 ± 0.001 

0.5 - 1.0 0.198 ± 0.002 ± 0.001 0.173 ± 0.006 ± 0.001 0.088 ± 0.026 ± 0.001 

1.0 - 1.5 0.209 ± 0.003 ± 0.001 0.191 ± 0.005 ± 0.001 0.095 ± 0.029 ± 0.001 

1.5 - 2.0 0.224 ± 0.003 ± 0.001 0.207 ± 0.006 ± 0.001 0.099 ± 0.036 ± 0.001 

2.0 - 2.5 0.233 ± 0.006 ± 0.002 0.215 ± 0.008 ± 0.001 0.100 ± 0.042 ± 0.003 

2.5 - 3.0 0.245 ± 0.008 ± 0.005 0.229 ± 0.011 ± 0.001 0.123 ± 0.052 ± 0.011 

3.0 - 3.5 0.264 ± 0.011 ± 0.009 0.228 ± 0.015 ± 0.001 0.080 ± 0.062 ± 0.017 

3.5 - 4.0 0.254 ± 0.015 ± 0.015 0.235 ± 0.020 ± 0.001 0.163 ± 0.080 ± 0.019 

Table 2 

Fully corrected production asymmetries of A, E, and 0, integrated over three x F re- 
gions. Errors are obtained by adding the statistical and systematic errors in quadra- 
ture. Results of ACCMOR [7,8] integrated over the range < x F < 0.35 are shown 
for comparison. 

E791 E791 E791 ACCMOR 

-0.12 < x F < 0.12 -0.12 < x F < < x F < 0.12 < x F < 0.35 

i, n t ii 0.207 ±0.001 0.242 ±0.002 0.124 ±0.002 0.119 ± 0.009 [71 

A — A L J 

A s -_ s+ 0.176 ±0.004 0.186 ±0.004 0.150 ±0.007 0.130 ± 0.050 [8] 

^n--n+ 0.099 ±0.013 0.098 ±0.018 0.100 ±0.021 0.107 ± 0.070 [8] 
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